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The time-evolution operator for the kicked Harper model is reduced to block matrix form when the 
effective Planck's constant h = 2tt^ and M and are integers. Each block matrix is spanned by an 
orthonormal set of A'^ kq (quasi-position/quasi-momentum) functions. This implies that the system's 
eigenfunctions or stationary states are necessarily discrete and periodic. The reduction allows, for 
the first time, an examination of the 2-dimensional structure of the system's quasi-energy spectrum 
and the study of, with unprecedented accuracy, the system's stationary states. 
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INTRODUCTION 



The Harper model was first introduced in 1955 to ap- 
proximate the dynamics of electrons confined to a 2- 
dimensional lattice while under the influence of a per- 
pendicular magnetic field 0. The model has since been 
generalised with the addition of a periodic kicking term 
so that it is classically chaotic under certain conditions 
0, Q ■ The dynamics of this kicked Harper model can be 
related to those of the kicked harmonic oscillator when 
a ratio of 1/4 exists between the kicking and oscillation 
frequencies. 

The kicked harmonic oscillator was originally proposed 
as a 2-dimensional model of charges moving in a homo- 
geneous static magnetic field while under the influence of 
an orthogonal time-dependent electric field 0, 0, IE 13 ■ 
The model is fundamentally different from other widely 
studied kicked systems because the natural frequency of 
the unperturbed system does not depend on energy. It 
therefore cannot be described using the KAM theorem 
. The quantized system has been proposed as a model 
for electronic transport in semiconductor super-lattices 
[lo| and for atom optic modeling in ion-traps |^ . 

The quantum kicked harmonic oscillator and the quan- 
tum kicked Harper model, in the same way as the quan- 
tum kicked rotator j ih^-Y be simulated efficiently on a 
quantum information processor [isl , using the Quan- 
tum Fourier Transform |l5| and Quantum Fractional 
Fourier Transform 16] . 

In this article we show how the evolution operator for 
Kicked Harper model may be reduced to block matrix 
form in the kq representation. This is then used to ex- 
amine the 2-dimensional structure of the systems quasi- 
energy bands as well as to construct the eigenfunctions 
of the evolution operator in the position and momentum 
representations. We also compare the structure of se- 
lected eigenfunctions with the classical phase-space using 
the Husimi quasi-probability distribution, see Fig. ^ 
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FIG. 1; A section of the Husimi distribution of a quantum 
stationary state of the kicked Harper model when K — L — 
1.5 and h — 2n/20. The distribution has a structure similar 
to the Poincare surface of section illustrated in Fig. |5| 



THE KICKED HARPER MODEL 



The generalised kicked Harper Hamiltonian is written 



as 



H = Wip) + V{q) J2 



(1) 



] = -oc 



It is usually assumed that W{p) — Lcosp and V{q) — 
K cos q although our analysis will apply to all functions 
that are 2tt periodic. 

The classical model displays some unusual phase space 
properties. In particular is the ordering of the Poincare 
surface of section into a lattice of stable elliptical cells 
when the perturbation is applied. Each cell is separated 
by a mesh of unstable trajectories which are collectively 
called the stochastic web, see Fig. |21 The web extends 
to all areas of the phase plane and grows in thickness as 
the perturbation strengths are increased 

Another interesting feature of the model is the exis- 
tence, under certain resonant conditions, of accelerator 
modes or islands. These modes allow for stable ballistic 
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FIG. 2: Poincare surface of section for the kicked Harper 
model when K = L = 1.5. 



motion throughout the phase space and are the main rea- 
son behind what is called anomalous diffusion or anoma- 
lous energy growth [sI IitIIisI . 

The quantum mechanical study of this problem has 
been dealt with in many guises and by many authors 
[IlllllIllllllllllllllMllilSllillS The delta 
function in means that the quantum mechanical time- 
evolution or Floquet operator may be factored as 



U{q,p) = P X Q = exp ( -\w{p)^ exp [~\v{q) 

(2) 

where p and q are the momentum and position operators 
whose eigenvalue equations are p\ p') = p'\ p') and q\ q') = 
q'\ q') respectively. The Fourier transform operation, 



e^p(^-q'p']f{q')dq', 



(3) 



provides the means to write ^ in either position or 
momentum representations. If it is understood that 
acts on everything to the right we may say that 
P{p) = T^P[q)T where P is a function of the operators 
p and q. We can therefore write the evolution operator 
U as 



U{q) = exp ( ~-W{q) ] J^exp 



-V{q) 



(4) 



The system's link with the Kicked Harmonic Oscillator 
can be seen from the fact that the the Floquet operator 
of that system can be written as 



UKHM{q,0) = exp 



-V{q) 



(5) 



where J--g represents what is called the fractional Fourier 
transform and is, up to a phase, the evolution operator of 



the simple harmonic oscillator in the position representa- 
tion |l3l l3l| . The ordinary Fourier transform is defined 
as !F = and time-evolves a state forward 1/4 of one 
full oscillation. 



THE FLOQUET OPERATOR IN THE kq 
REPRESENTATION 



The operators U and UKHM^'i-' f ) can both be shown 
to be simultaneously translationally invariant in posi- 
tion and momentum if /i = l-KMjN [H H^]. This 
suggests that the orthonormal kq (quasi-position/quasi- 
momentum) functions of Zak j32l |. which is also transla- 
tionally invariant in position and momentum, would be 
a convenient representation to work in. 

To define this representation we first introduce the 
Dirac delta comb or train as 



((7'|IIIi)^III,(q')^ E %'-"0- 



(6) 



n— — oo 



Then, in a similar way to Dana [27|, we consider the 
(slightly modified) kq functions, which we define as 



(g'|r, s, I) EE exp ^(g' - -) III,(g' - r) (7) 



with C„ 



In our notation the k of Zak is re- 



1 Y V27rS' 

placed by r and q is replaced by s so we may use q and 
q' as the position operator and variable without confu- 
sion. It has been shown that these functions form an 
orthonormal set with the labels 0<r<Z,0<s<a and 
a = 2-Khll m. 

In the momentum representation the states | r, s, Z) 
may be written as 



(p'|r,s,0 = .F((7'|r,s,0, 

1 /■ I I 



(8) 



2'nhJ \ ' ^'^'1^'^'^'*^^'' 

Cpexp('^(p'-^))lII„(p'-s), 



where Cp — 



/2-Rh 



It should be noted that the 



modification to Zak's definition is solely to preserve the 
symmetry between momentum and position representa- 
tions. 

The procedure for reducing expression to block 
matrix form is a technical exercise in labeling and re- 
ordering and therefore given in the appendix. There 
it is shown that for block reduction to be possible one 



3 



must choose the inversely related / and a so that they 
are both rational multiples of the periodicity of W and 
V. For h = 2ttM/N one may do this by either setting 
{I, a) = {2TrM/Ns,2Tr/Nr) or {I, a) = {2Tr/Ns,2TrM /Nr) 
with the added constraint that NrNg — N. The resulting 
block matrices are of dimension N x N and are spanned 
by the kets 



must be of the form 



I r + jh, s + kh, I) 



(9) 



where j € 0, 1, ...,Nr — 1 and A: e 0, 1, ...,-/Vs — 1- Impor- 
tantly, this set of kets defines the subspaces on which the 
system dynamics take place. 

As an example suppose H = Tr/2. One can make the 
choices {I, a) = (tt, tt), (27r, 7r/2) or (7r/2,27r) among oth- 
ers. The first choice (/, a) — (tt, tt) means the spanning 
subspaces are 



I r, s, tt) 

TT TT 

\r + ^,s+ -,tt) 



(10) 



for < r < ?i and < s < h. The second and third 
choices mean the subspaces are 



r, s, 27r) 
k + |,s,27r) 



and 



\r -\- Tr,s, 27r) 
\r + —,s,27r), 



l^'^+2'2> 



\r,s + Tr,-) 
I^'^ + T'2> 



(11 



(12) 



respectively, where again < r < h and < s < H. It 
can be shown that these three subspaces are actually one 
and the same, see 01 • In. the appendix, the elements of 
each block matrix are calculated and written in a concise 
form using the generalised Fourier matrices, see equations 
(|35|l and H37|) . As each block matrix is spanned by the 
kq states of solutions to the eigenvalue equation 



C/|V) = e'"|V'), 



^ J2 J2 i:^^k\r + jh,s + kh,l), (14) 



j=0 k=0 



where the ipj^k are the N elements of one eigenvector of 
either (|35|l or H37|) . This implies that generic eigenfunc- 
tions of U, with h — 2nAI/N, are translationally invari- 
ant in position and momentum. This is true regardless 
of the values of K and L. 

This result has important consequences for the inter- 
pretation of the systems dynamical behaviour, in particu- 
lar for the behaviour of initially localised states over long 
time scales. To measure the spread of a wavefunction 
over time one often examines the harmonic oscillator en- 
ergy function E = ^{q^+p"^)- One can classify the type of 
dynamics observed in terms of the growth rate of the en- 
ergy expectation value \E\ </>„) where | 0„) = (f)o)- 
A significant amount of effort has been made in re- 
cent decades to try and understand these growth rates in 
terms of system's quasi-energy distribution and the type 
of stationary states observed. The result above is impor- 
tant because it says that all eigenstates (at least when 
h = 2ttM/N) are extended and therefore suggests that 
qualitative differences in observed dynamical behaviour 
should depend on the quasi-energy spectrum alone. How- 
ever, this is not the full story. While the result deprives 
us the use of the eigenfunctions to explain the dynamics 
it provides us with another means, namely the size of the 
dynamical subspaces N. 

Analytical and numerical arguments by Borgovoni et 
al. '.2b\ have shown that iih = 27rp, where p is some irra- 
tional number, then the energy growth rate after a long 
time is linear. However they also show if h = 2ttM/N ^ 
the energy growth rate eventually becomes quadratic in 
time. This quadratic growth is sometimes known as 
quantum resonance and has been confirmed by a number 
of different numerical techniques [s^ . 

Although we cannot use the kq functions to examine 
the non-resonant situation h = 2'Kp (there being no way 
to choose {I, a) so that they both are rational multiples 
of 27r), we can always say that there exists a rational 
approximate M /N that is arbitrarily close to p and for 
which kq block reduction can be achieved. With this 
idea in mind it may be possible to relate the time before 
the onset of quadratic energy growth to the sizes N of 
the extended but finite subspaces. The constant linear 
energy growth rate observed for h = 2-Kp could then be 
explained as a direct consequence of the motion taking 
place on an extended but infinite space. 

The form of the Floquet eigenfunctions, H14|l . strongly 
suggests that the dynamical localisation observed in the 
kicked rotator j^^l does not occur in this model. Dynam- 
ical localisation is when the quantum expectation value 
(13) {4>n \E\(l)n) saturates around some constant value while 
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the energy of a classical ensemble, with the same initial 
marginals, continues to grow js^. The most widely cited 
explanation of this process is given by Gremple et al. 
[soi l . Here the phenomena is linked to Anderson localisa- 
tion ^3 by way of a mapping to a tight-binding system. 
However, a key point in this explanation is that the Flo- 
quet eigenfunctions are exponentially localised. Clearly 
the eigenfunctions of the kicked Harper model do not 
display this property. 

Indeed, the form of the eigenfunctions (|14|l suggests 
that, for all non-zero values of K and L, indefinite local- 
isation is unlikely. Note that after each operation of U, 
and because of differing quasi-energies, the interference 
between Floquet eigenfunctions that initially allowed the 
creation of a localised state, will gradually disappear. 
The fact that the dynamics take place on extended sub- 
spaces means that the wavefunction must eventually be- 
gin to spread. This should happen even if the initial state 
is placed inside a stable classical phase-space structure. 

The only way this dispersion or diffusion may be pre- 
vented is if the quasi-energies of the states that make up 
the superposition are all the same or only differ from each 
other by some rational multiple of 27r. It is difficult to 
envisage a situation where this could be arranged 




(a) uj± (r, s) with L = n, K = n 
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APPLICATION 



(b) u!^ (r, s) with L = n, K = n 



In this section some analytical and numerical exam- 
ples of eigensolutions for specific L and h are given. 
These are primarily included to instruct and to visually 
aid understanding. However they are also to demonstrate 
the consistency and validity of the results of the previous 
section and of the analysis given in the appendix. 



Explicit calculation with h = tt 

One can apply H35|l or (|37|l to the problem with H = n 
(M = 1 and N = 2). This has already been done in fs^l, 
however we may now give the results in terms of the two 
parameters r, s instead of just one. Using l|35|l it can be 
seen that 



U{r,s) = X{l,2)V{r,s)X{l,2)Q{r,s) (15) 

exp(~iL) cos{K) —i exp(—iL) sm{K)d \ 
-'iexp{iL)ain{K)9* exp{iL) cos{K) J 

where L = ^ cos s, K = j- cos r and — exp(ir/2). This 
calculation takes place on the subspace 




(c) aj±(r, s) with K = IQ,L - 
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k, j ' 5 -f /i, z) = 

but one could have also defined the subspaces 



(16) 



(d) uj+{r, s) with K = W,L = 11 



FIG. 3: u}±{r,s) ( see eq. HSU ) for h = tt and various values 
of K and L. 
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|r,s,0=(^oj' \r + f^^s,l)=i^^j, (17) 

and used equation (|37|l to perform a similar calculation. 
The eigenvalues of the matrix (|15(l can be seen to be 



e'^^i = cos(L) cos(if) ±iyi- cos^{L) cos^{K). (18) 

The surfaces uj±{r,s) for {K,L) — (tt, tt) and {K,L) — 
(10, 11) are plotted in Fig. |3| To emphasise that there 
may exist many different subspaces, labeled by r and s, 
that give rise to eigenfunctions with same quasi-energy co 
we have also illustrated some of the contour lines of cu^. 
The eigenvectors of ()15f) may be written as 

I ,+^ ^ ivr f _ exp(-iZ) sinO^)^ \ 

' ^ ' sin(L) cos(K) T Vl-cos2(L) cos^{K) J ' 

(19) 

where N± are normalizing coefficients. In the position 
representation these may be written as 



routines. As an example of this technique and the in- 
sight it provides into the system's quasi-energy structure 
we plot six quasi-energy bands as a function of r and s, 
for = i = 2 and S = 27r/400, in Fig. The figure 
demonstrates that the quasi-energy structure of the sys- 
tem, under these parameters, contains intersecting point 
spectra and continuous bands. 




FIG. 5: Some quasi-energy surfaces with K — 2, L = 2 and 
h — 27r/400. The spectrum contains intersecting discrete and 
continuous bands. 



=^o±o(g'|r,s,^)-K7Ao±i((7'|r,s + ft,7r), (20) 

where ip^f, is given by {r+jh,s + kh,TT\ip^) and {q'\r,s,l) 
is given by Q. An example of a single eigenvector, 
viewed from the position representation, is illustrated in 
Fig. H 

0.2 1 1 1 1 1 0.2 I 1 1 1 1 



o 
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q q 

FIG. 4: An eigenstate viewed from the position representa- 
tion, see II2UII and note '4T|. The stationary state is a super- 
position of two orthogonal kq functions. 



Convergents of ft / 2-kM/N 

Block reduction can only take place for h = 2ttM/N 
as there is no way to choose {l,a) so that they both are 
rational multiples of 27r when h = 2'Kp and p is irrational. 
However, it is possible to examine eigensolutions of the 
problem for rational approximates M/N of p. 

4 1 i 1 i 1 i 1 1 

2 




-2 - 
-3 - 



Numerical analysis with h = 27r/400 

In a similar way we may calculate the matrix elements 
of or |jS3) for general h = 2TTM/N. Explicit eigen- 
solution calculations for large N become unmanageable 
but can be done numerically using standard LAPACK 



-1l -10 -9 -8 -7 -6 -5 -4 

ln|pg - M/N| 

FIG. 6: The quasi-energy band structure with K = L = 1.5 
and h = 2-kM/N with M/N = 1/7, 1/8, 2/15, 3/23, 5/38, 8/61 
and 13/99 These are the first 7 convergents of the continued 
fraction expansion of pg — 1/ {6 + ag) where ag is the golden 
mean {V5 + l)/2. 
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Following |25| we examine convergents of the continued 
fractions expansion of h = 2-Kpg where pg = 1/(6 + cjg) 
and (Tg is the golden mean (-s/S + l)/2. The results are 
given graphically in Fig. The denominator N of each 
convergent increases as we approach p and therefore so 
must the number of quasi-energy bands. Recent analysis 
suggests that as M/N — > p „ the sp ectral distribution may 
even become multi- fractal |2J, |22| . 

As an example of a stationary state of a rational ap- 
proximate to irrational pg we choose M/N = 233/1775. 
We plot one of the calculated eigenstates in the position 
representation in Fig. [T] At a glance the eigenfunction 
first appears to be continuous but on closer inspection 
it can be seen to be discrete. Importantly, this eigen- 
state and indeed all eigenstates calculated for a rational 
approximate to irrational p are still periodic and still ex- 
tended. There is no need to give them extra significance. 
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(a) Phase space plot of the 
accelerator modes responsible for 
anomalous diffusion with 
K = L = 2iT + .1 
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0.1 0.2 



FIG. 7: An eigenfunction in the position basis with K = 1 and 
L — 7 with h = To generate this state a superposition 

of the A'" — 1775 kq states {q'\0, 0+jh, 2tt) with j running from 
to A'' was used, note 



-2 2 4 6 



(b) The Husimi distribution of a 
stationary state for 
K = L = 2iT + O.l and h = 27r/200 



Quasi-probability distributions 

As already mentioned one may reconstruct each state 
in the position or momentum basis by using equations 
(0) or 0. However, we may also examine the Wigner 
and Husimi distributions of these states. As a brief ex- 
ample of this technique we plotted in Fig. ^the Husimi 
quasi-probability function of a stationary state that has 
the same web-structure as the classical Poincare surface 
of section illustrated in Fig. |21 The translational invari- 
ance of this stationary state is clearly visible although 
the Husimi distribution obscures it's discretised nature. 

Now, we also briefly examine a stationary state re- 
sponsible for quantum anomalous difi^usion • Classical 
anomalous diffusion is due to stable phase space struc- 
tures occurring for particular values of K and L. When 
K ~ L the stability conditions for the existence of these 
structures may be written as 2mT < K < y^2ri7rp + 4 
for integer n [Sg- The plot in Fig. [Sjis evidence of a sort 



FIG. 8: Accelerator modes. The Husimi distribution of this 
stationary state has peaks pri marily around the classical ac- 
celerator modes. Note that Fig. |8(a)| is not a Poincare surface 
of section but rather a plot of the initial conditions which, 
after 800 iterations of the classical mapping, have energies 
> 3500 



of scaring as classical points initially on an accelerator 
mode never return to the same region. 

These examples illustrate how the quantum eigensolu- 
tions may resemble classical phase space structures. In 
further studies this may allow one to clearly distinguish 
between quasi-energy distributions associated with sta- 
ble classical trajectories and those related to the unstable 
trajectories. 
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CONCLUSION 

The Floquet operator for the kicked Harper model has 
been redueed to block matrices for h — 2'kN /M and pos- 
itive integers M and N . As each block is spanned by a 
set of N orthogonal kq (rs) functions, any eigenfunction 
of the Floquet operator must be a superposition of this 
set and is therefore periodic in position and momentum. 

We have explicitly calculated the block matrix ele- 
ments and showed that the matrices themselves may be 
written concisely using the generalised discrete Fourier 
transform. Using this method we have performed anal- 
ysis on the system for various values of the quantum 
parameter h and the parameters K and L. We have 
illustrated, for the first time, the 2-dimensional struc- 
ture in the system's quasi-energy spectrum and demon- 
strated, albeit with just one example, that the distribu- 
tion may be made up of intersecting discrete and contin- 
uous spectral bands. The behaviour of the spectrum as 
h approaches an irrational multiple of 27r is also demon- 
strated. 

In addition to this, selected stationary states of the sys- 
tem have been calculated. It is demonstrated that they 
are, by construction, periodic and discrete in momentum 
and position. Two examples of stationary states whose 
Husimi distributions resemble classical phase space struc- 
tures have also been provided. These results should help 
open up the way for a more detailed study of the system's 
eigensolutions and their relationship with the classical 
system. 

ACKNOWLEDGMENTS 

The author sincerely thanks D. M. Heffernan for the 
many fruitful discussions relating to this work. 

APPENDIX 

The elements of the Floquet matrix in the kq repre- 
sentation may be written as 

(r",s",/|C/|r',s',0= ^(r",s",;|P|r,s,0 (21) 

r.s 

x(r,s,/|Q|r',5',0 (22) 
We first calculate the expansion 

Qr,sy,s' = {r,s,l\Q\r',.3',l) (23) 
= J J {r,s,l\q'){q' \Q{q)\q"){q"\r\s',l)dq'dq". 

By substituting (0 into (|23|l and after some manipulation 
it can be seen that non-zero matrix elements occur only 
if r = r'. We therefore write 



Qr,s,r,s' = C'^^exp f-^(^ + nO j X (24) 
^Q(r + mZ)expr^(^+mOj . 

It is now that the choice of {I, a) becomes important for 
the block reduction. Setting {I, a) = {2i:M /Ns,2t: /Nr) 
one can see after some consideration, that the summa- 
tions are only well-defined if s — s' = (j — k)h with 
j, k = O...Ns — 1. To this end we make the substitutions 
s = s + kh and s' = s + jh to get 

Qr.s = C^|^exp(-i(fc+^)(J + n)?) X (25) 

n 

^Q{r + ml) exp + ^){^^ + m)l^ . 

m 

A simple analysis shows that the two infinite summations 
can be reduced to finite ones with integers n and m run- 
ning from to A'^s — 1. By representing the orthonormal 
kq kets as 

\r,sj) = (1,0,0,0,. ..,0,0)^ 

\r,s + h,l) = (0,1,0,0, ...,0,0)^ 
\r,s + 2h,l) = (0,0, 1,0,. ..,0,0)^ 

|r,s-f (A^,-l)ft,0 = (0,0,0,0,...,0,l)^, 

one can write these finite summations as matrix multi- 
plications. That is 

«('-'"='^.^^(i'^)^(''>''G'5)' <^'" 

where D is the diagonal matrix defined by 

I?(n-,=exp(^--iy(r + j7)^, (27) 
and F(a,b) is the generalised Fourier matrix 

F,- = c^(^+'^)('=+^\ (28) 

with ujq = exp{-i2nM/Ns)/^/^^. 

A similar analysis can be done for the matrix elements 
(r", s", I \P\ r, s, I) by expanding in the momentum basis 
giving 
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where in this case the matrix D is defined as 

D{W),,j = exp (-'j^W{s + ja)^ , (30) 

and Up = exp{—i2TT/Nr)/^/Nr. Of course to write it Hke 
this the kets spanning the subspaces are represented as 



where of course j x k = N. If we retain the basis labeUng 
of (|31|) we must rearrange P{s,r) to get 

V'ir, s) = X{Nr, Ns)^rir, s)X{Nr, N,). (34) 

The total evolution operator, (0), on the r,s subspace, 
can then be written as: 



\r,s,l) =(1,0,0,0,...,0,0)^ 
\r + h,s,l) = (0,1, 0,0,.. .,0,0)^ 
\r + 2h,s,l) = (0,0,1,0,. ..,0,0)^ 



U{r,s) =r'ir,s)Q{r,s) 



(35) 



However, we could just as easily have retained the basis 
labeling of (p?^ by reordering Q(r, s) to get 



|r + (7V^-l)S,s,0-(0,0,0,0,...,0,l)^. 

It should be noted here that the expression for uip does 
not contain the factor M. However, this is only a 
consequence of our choice {l,a). If one sets {I, a) = 
{2n/Ns,2'KM/Nr) we get Ug = exp{~i2n/Ns)/ y/Ws and 
ujp = exp{-i2TTM/Nr)/y/N^. 

Together Q{r, s) and F(r, s) act only on the subspaces 
spanned by the kets | r + jTi, s + fcfi,, I) for j = Q...Nr — 1 
and k = O....Ns — 1. We define the the block diagonal 
matrix Q as 



Q(ri,s) 
g(r2,s) 














Qi'rN^,s) 



,(31) 



where Vn — r + (n — l)h and each block matrix Q{rn, s) 
is given by (|26|l . This matrix calculates the operation 
Q{q) = eTi"^('' over the whole subspace. In the same 
way we define the block diagonal matrix V as 



V(r,s) 



P{r,si) 
P(r,S2) 
















,(32) 



for the e"^'^-* operation. Here, s„ = s + (n — l)h and 
each block matrix P{r, s„) is given by (|29|) . It is helpful 
to note here that C^Cp — 1. 

In order to multiply these block matrices together one 
must first re-order the elements of one of the matrices. 
This can be accomplished through the use of a simple 
swap matrix. If n is the n*'' column vector of the N x N 
identity matrix we can define the swap matrix X{j, k) as 



X(i,fc)==[l,l+i, 



l + (fc-l)j, 2, 2+j, ... 
j + (fc-2)i,j X fc], (33) 



Q'(r, s) - X{Nr, Ns)Q{r, s)X{Nr, 7V,)t, (36) 
and now the Floquet operator is written as 



Uir,s)=rir, s)Q\r, s). 



(37) 



It is the special case of when Ng = 1 and Nr ^ N 
that gives us method used in jl^ I^Ol to analyse this 
system. The matrix Q{r, s) reduces to one with just 
exp(— jy(r„)/?i) along the diagonal and V{r,s) is given 
by H29|) . In total we have 



f/(r,s) = V{r,s)Q'{r,s) =V{r,s)Q{r,s) , 

:W{s,,)) (38) 



'2a h' 



Equation (|39|l bears a striking similarity to the position 
representation of ?7 in Q), albeit with discrete values of q. 
Understandably, this often led to the interpretation that 
the eigenvector elements of U{r, s) are discrete values of 
the actaa/ eigenfunctions of U in the position basis. How- 
ever, as this analysis has just shown, each block matrix 
is spanned by the vectors 



\r + jh,s + kh, I) 



(39) 



where j £ 0, 1, ...,Nr - 1 and fc e 0, 1, Ns - 1. There- 
fore solutions to the eigenvalue equation, U\ t/') = e'"| 
must be of the form 



IV) = ^ 5] V',,fek + j^,s + fcft,0- (40) 

j=0 k=0 

An eigenstate viewed from the position representation 
{q'\ip) is therefore expanded as 
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(9'I^)=E E + (41) 

j=0 k=0 

where -i/jj.fc — {r + jh, s + kh,l\'ip) and {q'\r + jh, s + kh,l) 
is defined bylQ. By construction the values ipj^k are the 
individual elements of the eigenvectors of the matrices 
C/(r,s). 
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